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Параметрический синтез системы управления 
электроприводом манипулятора с учетом 
неопределенности момента инерции нагрузки 


В статье рассматривается корневой подход к расчету параметров системы автоматического управления 
звеном манипулятора промышленного робота при условии значительных изменений момента инерции 
нагрузки, линейно входящего в коэффициенты характеристического уравнения системы. Задача 
решается путем размещения семейства корней системы в заданной области О, называемой областью 
качества. В результате на основе применения полей корневых траекторий определяется радиус 
дисковой области значений неопределенного параметра, обеспечивающих О-устойчивость системы. 


Введение 


При функционировании сложных технических объектов, как правило, наблюдаются 
отклонения параметров от расчетных значений, которые могут быть значительными, 
и по этой причине могут оказывать отрицательное воздействие на качество и устойчи- 
вость работы объектов. Следовательно, большое значение имеет разработка эффективных 
методов построения систем автоматического управления (САУ) оборудованием и, в 
частности, приводами манипуляционных роботов, которые позволили бы гарантировать 
требуемое качество работы указанных устройств в условиях неопределенности [1-3]. 
Решить рассматриваемую задачу можно с помощью самонастраивающихся систем [4] 
и методов обеспечения робастности. Методы построения робастных систем, как правило, 
базируются на использовании алгебраического [5], [6], частотного [7] и корнево- 
го [1], [8] подходов. 

Целью данной работы является применение корневого метода параметрического 
синтеза динамических систем, построенного на основе использования полей корневых 
траекторий кругового образа [1], [8] при синтезе системы управления приводами 
манипуляционного робота с учетом изменяющегося момента инерции нагрузки. Указан- 
ные системы позволяют обеспечить требуемое качество управления путем размещения 
корней характеристического уравнения в заданной области. 


Постановка задачи 


Рассмотрим систему подчиненного управления приводом одной степени под- 
вижности манипулятора, структурная схема которой приведена на рис. 1. Полагается, 
что схват манипулятора совершает относительно медленные движения, при которых 
эффектами взаимовлияния между всеми его степенями подвижности можно пренебречь. 
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Рисунок 1 — Структура системы управления положением звена манипулятора 


На рис. 1 И’.(5) = И/.(5) 5, где И’($) — передаточная функция объекта управления; 
Кс и К» - коэффициенты усиления регуляторов по скорости и положению; Т - по- 
стоянная времени ПИ-регулятора; (И, — задающее напряжение; ф — угол поворота 
вала привода; ф. — желаемое значение угла поворота; © — угловая скорость вращения 
выходного вала привода. 

Передаточная функция объекта управления имеет следующий вид: 


о (1) 


й О АВЕ з ре В: 29 
ь +1) — 5 +в +1) — 5 +С,5 
(ето С (ов + Ли) С ь 


м м 


где К,, Г,- активное и индуктивное сопротивления якорной цепи электродвигателя 
соответственно; /, — приведенный момент инерции нагрузки; /», — момент инерции 
якоря электродвигателя; Си Суи- конструктивные коэффициенты электродвигателя; 


К и К, -— коэффициенты усиления обратных связей по скорости и по положению 
соответственно; 5 — символ преобразования Лапласа. 

С учетом выражения (1) характеристическое уравнение системы управления 
электроприводом манипулятора имеет вид: 


4; а 53+ С Си м СиК!К. 5+ СиК›К,К. = 
Г (ов + ЛГ о + ЛОТ О + ЕТ 

В зависимости от массы переносимых изделий и текущей конфигурации мани- 
пулятора величина /, может изменяться в широких пределах [4], [9], что приводит к 
изменению коэффициентов характеристического уравнения (2). 

В плоскости корней 5 = в + 1» уравнения (2) зададим область О в форме трапеции, 
называемую областью качества. Эта область ограничивается двумя линиями равной 
степени устойчивости: Г.’ с уравнением о = п’ и Г.” с уравнением с =1"'; а также 
двумя линиями постоянного демпфирования: [+в и [вс уравнениями © = + ^о, где 
^ = & +В (В - угол наклона линий [+в и Г вкоси 0), что равносильно заданию 
пределов изменения показателей качества системы — степени устойчивости 1 и 
колебательности В. Учитывая наличие существенной неопределенности параметра /», 
ставится задача нахождения области О значений /„, при которых неизменно обеспе- 
чивается расположение корней уравнения (2) внутри заданной области качества 0. 
Решение этой задачи позволит удержать качественные характеристики т и В в уста- 


(2) 
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новленных пределах, гарантируя тем самым О-устойчивость системы и выполнение 
условия 


ЕР 53 О, (3) 


где / = 1, п . Задача будет решаться с использованием полей корневых траекторий. 


Описание метода синтеза 


Опишем рассматриваемую динамическую систему характеристическим урав- 
нением вида 


ая... Ча я+а=0, (4) 


где а; — коэффициенты уравнения, / = 1, и. Рассмотрим случай, когда один из пара- 
метров системы (обозначим его через А), линейно входящий в коэффициенты ее 
характеристического уравнения, является неопределенным. Выделим этот параметр 
К путем преобразования уравнения (4) к виду функции комплексного переменного 
Е =} (5) = я и(о, 0) +1(0,0), (5) 
(5 
где $ф(5) и \у(5) — некоторые полиномы от комплексного переменного 5; и(о,65), у(6,6) — 
гармонические функции независимых действительных переменных с и 6. 
Выражение (5) представляет собой уравнение корневого годографа системы в 
общем виде, которое позволяет отображать некоторые образы (кривые) произволь- 
ной геометрической формы, заданные в плоскости неопределенного (свободного) 
параметра А, на плоскость комплексного переменного 5 (плоскость корней). Введем 
определения для ряда используемых далее понятий. 
Пусть ф(и,у) = С -— кривая в плоскости А, которую назовем образом корневого 
годографа. Тогда, используя соотношение (5), получим уравнение 


Ес, 0) = С, (6) 


где Е(о,6) = ф(и(о,6),у(с,0)), С — некоторая постоянная для конкретного образа ве- 
личина, которую назовем параметром образа корневого годографа. Тогда функция 
Е= Е ,6) (7) 

представляет собой функцию поля корневых траекторий [1], [10]. Очевидно, что урав- 
нение (6) является уравнением линий уровня поля корневых траекторий. 

Определение 1. Полем корневых траекторий динамической системы с харак- 
теристическим уравнением (4) назовем совокупность кривых (6), когда параметр С 
образа корневого годографа изменяется в промежутке (0, - оо). 

Рассмотрим образ корневого годографа в форме окружности (круговой образ), и 
с целью формирования поля корневых траекторий кругового образа [1] запишем урав- 
нение корневого годографа кругового образа 


Е(о,6) = (и(с,6) — а)? + (у(о,6) — Б)! = С, (8) 
где С=р», р — радиус окружности-образа, заданного в плоскости свободного пара- 
метра К; а, Б — координаты центра окружности-образа по осям и и у соответственно. 
Уравнение (8) представляет собой уравнение линий уровня поля (7) для поля корне- 
вых траекторий кругового образа, а радиус р является параметром образа, который 
назовем также параметром поля в случае, если 0 <р<- <. Ввиду конформности 
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отображения, реализуемого функцией, обратной функции (5) [8], очевидно, что линии 
уровня поля кругового образа формируются в виде замкнутых кривых (годографов) 
в количестве для каждого годографа, не превышающем и, а центр окружности-образа 
поля отображается в виде точек в количестве и для каждого годографа, располагаю- 
щихся внутри этих кривых. Следовательно, изменяя расположение образа корневого 
годографа, можно задавать ориентацию отмеченных замкнутых кривых, размещая 
тем самым ограниченные ими области корней желаемым образом, в частности, в за- 
данных в плоскости 5 областях различной формы. 

Определение 2. Точки в плоскости 5 комплексного переменного, определяемые 
посредством отображения центра кругового образа, заданного в плоскости свобод- 
ного параметра, на плоскость 5, реализуемого с помощью функции, обратной к (5), 
назовем центрами локализации поля корневых траекторий кругового образа. 

Определение 3. Локальными назовем линии уровня поля корневых траекторий 
кругового образа, ограничивающие замкнутые односвязные области, отображаемые с 
помощью функции, обратной к (5), и не содержащие в качестве внутренних или гра- 
ничных точек точки ветвления этой функции. 

Определение 4. Глобальными назовем линии уровня кругового поля корневых 
траекторий, ограничивающие замкнутые односвязные области, отображаемые с помощью 
функции, обратной к (5), и содержащие в качестве внутренних или граничных одну 
или более точек ветвления этой функции. 

На основании отмеченных выше закономерностей определим искомую область 
р значений неопределенного параметра А в форме круга (диска) в комплексной 
плоскости А, ограниченного окружностью-образом радиуса р, и с помощью функ- 
ции (5) отобразим его таким образом, чтобы полученное при этом поле корневых 
траекторий (7) полностью располагалось в пределах заданной области качества О и 
было бы по возможности максимальным. С этой целью осуществляется вписывание 
линий уровня поля корневых траекторий кругового образа в область качества и 
определяется такая линия [, которая ограничивает область корней К, полностью 
принадлежащую области О и отображаемую с помощью (5) на плоскость А в форме 
диска-образа Р максимально возможного радиуса р = риах. Указанное вписывание 
выполняется посредством определения точек касания # (/= 1, т, где т — число точек 
касания) линий уровня поля и границ области 0. Очевидно, что линии уровня могут 
быть вписаны в область О только в случае расположения центров локализации поля 
внутри этой области, то есть когда это поле определенным образом ориентировано. 

Рассмотрим основные этапы решения задачи. 

1. Ориентация поля. Предположим, что центр окружности-образа находится на 
действительной оси и. В этом случае центры локализации поля могут быть размещены: 
— в нулях функции (5), т.е. в начальных точках корневого годографа Теодорчика — 

Эванса [1] исследуемой системы управления; 
— на участках положительных ветвей КГТЭ системы, расположенных в пределах 
области О. 
Координаты и = аи ь = В соответствующего центра образа определяются 
на основе формулы (8). 

2. Вписывание поля по горизонтали, то есть в область, ограниченную только 
заданными линиями равной степени устойчивости. При этом будет использован 
градиент поля корневых траекторий 


(9) 


эта4Е (с, 6) = ее ©) р о: 


{е [6.63] 
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где 1,/ — проекции единичного вектора, направленного по нормали к линии уровня 


поля. Проекции градиента на оси координат с и ® — эта4. Е(с,6) и ета4ь Е(с,0) соот- 
ветственно равны 


ОЕ(в, 6) 
до ` 


Поскольку в искомых точках касания проекции градиента на ось {равны нулю, 
определим координаты этих точек, составив две системы уравнений: 


эта4Е(о,6) = 


О. эта „Е(о,60)= 
до 


ОЕ (с, 6) ый 
00 ; (10) 
с=1' 
ОЕ(с, 6) 0 
00 (11) 
б=п" 


3. Вписывание поля по вертикали, то есть в область, ограниченную только за- 
данными линиями постоянного демпфирования. 

Координаты искомых точек касания определятся из следующей системы урав- 
нений: 


ОЕ(с, 6) 
06 


(-с) (-0)=0 


Ф=ЕЛО 


ОЕ(в, 6) 
а (12) 


где первое уравнение представляет собой уравнение касательной к линии уровня поля, 
проходящей через начало координат. 

С целью определения дисковой области О значений неопределенного парамет- 
ра А максимального радиуса риах для каждой точки касания, найденной путем решения 
систем (10) — (12) с использованием формулы (8), вычисляется значение радиуса 
образа р, и среди найденных таким образом значений этого радиуса выбирается ми- 
нимальное р = р, которое представляет собой радиус риах = р искомой области 2, 


удовлетворяющей условию (3). 


Обеспечение требуемых показателей качества САУ 


Для расчета значений неопределенного параметра, при которых будет обеспечи- 
ваться требуемое качество САУ, зададим значения постоянных параметров рассмат- 
риваемой системы и, подставив их в характеристическое уравнение (2), перепишем 
это уравнение в виде: 


$ +5505° + А(5’- 3605 + 28800) = 0. 

Поскольку величина /„ изменяется в широких пределах, то неопределенный (сво- 
бодный) параметр определим в виде А = И(/, +/,). 

Вводимая область качества ограничивается линиями равной степени устойчивости 
с =[' =- 65, в =[;’ =- 400 и линиями постоянного демпфирования [+в, для кото- 
рых Л, = (+В) = +07 (рис. 2). При этом используются полиномы $(5) =5' + 5505, 
ли(5) = 5’ + 3605 + 28800. В соответствии с приведенной выше методикой область О 
значений А ищется в форме круга (диска) радиуса р = риаах. Поскольку локализация поля 
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определяется конфигурацией годографа КГТЭ, выполним его построение и определим 
полюсы передаточной функции разомкнутой системы (нули функции (5)): ру =- 550, 
2 = рз = ра = 0 и ее нули (полюсы функции (5)): 21 =—120, 2 =—240 (рис. 2). 


И КСТЕ НТ = 123] 
Ее ЕЗ \Уем» Нер Ввод данных рт | 


Г Задание полиномов: Шир. рис. Выс. рис. 


С Коэффициенты Г 4 Е 


(® Корни 


<® Построение 


Г Данные: 
Сохранить Загрузить] С Вписывание 


пит: |500  зер СХ а: р 
пах; |130 з{ер и: |1 Ь: р 
чел: [55 — кой: [07 в [000 


де: |-400 Кор2: |-0.7 Г Область 


Г Диагностика м ктэ 


Отам Сапсе! рев аз зейез | 


г Числитель: 


Количество 


Р Г Мнимая часть 


№| Действительная 
1] 120 
2 


-240 


Г Знаменатель; 
Количество 


В Г Мнимая часте 


Рисунок 2 — Расположение корневого годографа Теодорчика — Эванса 
относительно области качества О 


Соответствующие полюсы на рис. 2 отмечены крестиками, а нули — кружками. 
Очевидно, что все полюсы выходят за пределы области О и, как отмечалось выше, в 
этом случае вписывание поля невозможно. По этой причине изменим конфигурацию 
КГТЭ путем «настройки» полюсов таким образом, чтобы они расположились в пределах 
данной области. Для этого зададим следующие значения полюсов: р! =-—100, р› =—140, 
рз=-170, р4 =-290 (рис. 3). При этом характеристическое уравнение рассматриваемой 
системы принимает вид: 


9+ 7005 + 173,7-10357' + 182,72-105 + 690,2-109+ &(5? + 3605 + 28800) = 0. 


В результате указанной «настройки» часть ветвей КГТЭ также расположилась в 
пределах области О и вписывание линии уровня стало возможным. Согласно ранее 
описанному методу, применим круговое поле корневых траекторий. С этой целью 
вначале выполним построение функции поля Е = Ё(с,6) (7) для заданной системы 
исходя из предположения, что центры локализации поля располагаются в полюсах р; 
(1 = 1,4) (эти полюсы обозначены крестиками), т.е. центр окружности-образа находит- 


ся в начале координат плоскости К. 
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ШИ КСТЕ | =[51 >| 
Е Е СИИ > 


Задание полиномов: Шир. рис. Выс. рис. 
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Рисунок 3 — Поле корневых траекторий кругового образа, вписанное 
в заданную область качества О 


Соответствующее поле на рис. 3 показано четырьмя линиями уровня: / (сфор- 
мирована локальными линиями /', /"" и глобальной линией /\"), 2 (сформирована 
локальной линией [' и глобальной линией [5"), [3 (сформирована локальной линией /з' 
и глобальной линией /[5") и [4 (глобальная линия). Линия /1 ограничивает трехлистную 
область корней, линия [5 — двулистную область корней КА, линия [3 — также двулистную 
область, а линия /4 — однолистную область корней. Далее на основе систем уравнений 
(10), (11) и (12) выполним вписывание указанных линий уровня в заданную область 0. 
В результате решения полученных систем уравнений определим значения радиусов 
р: ({= 1,4) дисковой области О значений параметра А: р: = 6443,18; р> = рз = 9607,52; 
р4 = 44049,11, соответствующие линиям р", [5"' и [1 на рис. 3 и соответственно точкам ка- 
сания № (радиус рт), В (радиус р> = рз) и & (радиус ра). Искомой является область, 
ограниченная окружностью, максимально возможный радиус которой р = рьах опреде- 
ляется минимальным р из полученных значений р; (1= 1,4): р = риах = пайа(ру, ро, рз, ра) = 


= р! = 6443,18. Эта область отображается на плоскость 5 в виде области корней К, 
ограниченной линией [5”, касающейся границы области О в точке Ь. Все значения (, 
лежащие внутри диска О радиуса р < 6443,18, обеспечивают требуемое расположе- 
ние корней в области О (3), а следовательно, и соответствующее качество процесса 
управления. 
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Заключение 


На примере расчета параметров системы автоматического управления электро- 
приводом одной степени подвижности многозвенного манипулятора, у которого в 
широких пределах изменяется приведенный момент инерции, линейно входящий в 
коэффициенты характеристического уравнения всей системы, показано эффективное 
применение корневого полевого метода размещения корней этой САУ в заданной 
области О. Посредством вписывания линий уровня кругового поля корневых траек- 
торий в трапецеидальную область качества О определяется радиус дисковой области 
значений неопределенного параметра системы, неизменно гарантирующих ее О-ус- 
ТОЙчивость. 
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А.А. Несенчук, В.Ф. Фларетов 

Параметричний синтез системи управлння електроприводом маншулятора з урахуванням 
невизначеност! моменту 1нерцИ навантаження 

У статт! розглядаеться кореневий шджд до розрахунку параметрв системи автоматичного управл!ння 
ланкою маншулятора промислового робота за умови значних зм1н моменту 1нерци навантаження, ли но 
вхдного в коефиценти характеристичного равняння системи. Задача розв’язуеться шляхом розмицення 
с1мейства кор!ння системи в задав област! О, званою областю якост!. В результат! на основ! 
застосування полв кореневих траектор!й визначаеться радлус дисковоГ област! значень невизначеного 
параметра, забезпечуючих О-спйюсть системи. 


АПа А. №езепсйий, Иаатиг Е. ЕЙагеоу 

Рагатейле Зупез$ оЁ Фе Сопи`о! Зу$ет ог Машрша юг Шесыле Омуе Сопз4етто пе Гоа 
МотепЕ оЁ шегНа Опсег{айциу 

ш Фе рарег Фе гооё 10си$ арргоасб 15 сопз14егеЯ №0 са]си]айп» рагатеегз оЁР Фе шдизла| гобой 
таприаюг сопёго] 1оор ш соп4! оп оЁ забуапна| уапайоп оЁ Фе 1оа4 шегна поглеп{ ещегте Ппеау 1е 
зу$ет свагасег1$Ис ефааНоп соеЁЯс1еп5. ТВе 1азК 15 зо]уе4 Ъу р!асетепе оЁ Фе сопго| зубет гос 
атПу улет Фе этуеп доташ О, Ус 15 патед Фе диаШу дотат. ЕтаПу, оп Фе Баз1$ оЁе гос [0си$ 
Ве!4$ аррПсайоп, Фе га@па$ оЁ {Ве 91с-5Варе4 гезлоп оЁ Фе ипсецатш рагатейег уа[аез 15 Чеегтлтеа \/БлсЬ 
эпагащее фе зубет О-5а у. 
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